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>; ABSTRACT 

S: 

We present the first case in which the BEER algorithm identified a hot Jupiter 
in the Kepler hght curve, and its reality was confirmed by orbital solutions based 
on follow-up spectroscopy. Kepler-76b was identified by the BEER algorithm, 
which detected the BEaming (sometimes called Doppler boosting) effect together 
with the Ellipsoidal and Refiection/emission modulations (BEER), at an orbital 



^' 

P4 period of 1.54 days, suggesting a planetary companion orbiting the 13.3 mag 

O ■ F star. Further investigation revealed that this star appeared in the Kepler 

-t— > ■ eclipsing binary catalog with estimated primary and secondary eclipse depths of 

c^ ■ 5 X 10"'^ and 1 x 10~^ respectively. Spectroscopic radial- velocity follow-up obser- 

vations with TRES and SOPHIE confirmed Kepler- 76b as a transiting 2.0 ±0.26 
^ ■ Mjup hot Jupiter. The mass of a transiting planet can be estimated from either 

the beaming or the ellipsoidal amplitude. The ellipsoidal-based mass estimate 
00 ■ of Kepler-76b is consistent with the spectroscopically measured mass while the 

^ ■ beaming-based estimate is significantly infiated. We explain this apparent dis- 

Q ■ crepancy as evidence for the superrotation phenomenon, which involves eastward 

cn ■ displacement of the hottest atmospheric spot of a tidally-locked planet by an 

equatorial super-rotating jet stream. This phenomenon was previously observed 
only for HD 189733b in the infrared. We show that a phase shift of 10.3 ±2.0 de- 
grees of the planet reflection/emission modulation, due to superrotation, explains 
the apparently infiated beaming modulation, resolving the ellipsoidal/beaming 
amplitude discrepancy. Kepler- 76b is one of very few confirmed planets in the 
Kepler light curves that show BEER modulations and the first to show superro- 
tation evidence in the Kepler band. Its discovery illustrates for the first time the 
ability of the BEER algorithm to detect short-period planets and brown dwarfs. 
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Subject headings: binaries: spectroscopic — methods: data analysis — planets 
and satellites: detection — stars: individual (Kepler-76, KIC 4570949) 



Introduction 



CoRoT and Kep ler have produced hundreds of thousan ds of nearly uninterrupted high 
precision light curves (JAuvergne et all2009l : lKoch et al.ll2010l ) that enable detection of minute 
astrophysical effects. One of these is the the beaming effect, sometimes called Doppler boost- 
ing, induced by stellar radial velocity. The effect causes a decreas e (increase) of the brightnes s 
of any light source receding from (approaching) the observer (JRybicki fc Lightmanlll979h . 
on the order of 4fr/c, where Vr is the radial velocity of the source, and c is the velocity 
of light. Thus, periodic variation of the stellar radial velocity d ue to an orbiting corn pan- 
ion produces a periodic beaming m odulation of the stellar flux. iLoeb fc Gaudil (120031 ) and 
Zucker. Mazeh &: Alexander! (120071 ) suggested to use this effect to identify non-eclipsing bi- 
naries and exoplanets in the light curves of CoRoT and Kepler. The precision of the two 
satellites is needed because even for short-period binaries, with large radial-velocity (RV) 
orbital amplitudes, the beaming effect is small, on the order of 100-500 ppm (parts per 
million). 

As pred icted, several studies identified the beaming effect in short-period known eclips- 



ing b i naries (Ivan Kerkwiik et al. 



2011: Bloemen et al. 



2011 



2012 



2010 



Breton et al. 



2012 



Rowe et al.ll201ll:ICarter et al.ll201ll : lKipping fc Spiegel 



; I Weiss et al.l 



20121 ). Yet, space missions 



data can be used to i dentify non- eclipsing hmaiies thro ugh detection of the beaming effect 
( JLoeb fc Gaudill2003l : IZucker. Mazeh fc Alexandeiil2007l ). However, the beaming modulation 
by itself might not be enough to identify a binary star, as per iodic modulations could be 
produ ced by other effects, stellar variability in particular (e.g 
2004h . 



Aigrain. Favata &: Gilmore 



To overcome this problem, the BEER algorithm (JFaigler fc Mazehll201ll ) searches for 
stars that show in their light curves a combination of the BEaming effect with two other 
effects that are produced by a short-period co mpanion — t he Ellipsoidal and the Reflection 



modulations. The ellipsoidal variation (e.g., iMorris 



between the two components (see a review by iMazehl (120081 )). while the reflection/heating 



1985) is due to the tidal interaction 



variation (referred to herein as the reflection modulation) is caus ed by the lurninosity of each 



component that falls on the facing half of its companion (e .g., IWilsorull990l : iMaxted et al. 
2OO2I : iHarrison et aPboOsi IPor et al.ll201ol : iReed et al.lboiol ). Detecting the beaming effect 
together with the ellipsoidal and reflection periodic variations, with the expected relative 
amplitudes and phases, can indicate the presence of a small non-eclipsing companion. Re- 
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cently iFaigler et al.l ( 120121 ) reported RV confirmation of seven new non-ecUpsing short-period 
binary systems in the Kepler field, with companion minimum masses in the range 0.07-0.4 
Mq, that were discovered by the BEER algorithm. 

For brown-dwarfs or planetary companions the beaming effect is even smaller, on the 
order of 2-50 ppm. Interestingly, several studies were able to detect this minute effect in sys 



tems with transiting brown dwarfs and vlanets (IMazeh fc Faiglerll2010l: 



Mazeh et al.ll2012l : I Jackson et al.ll2012l : iMislis et al.l 120121 : iBarclav et al. 



Shporer et al.l 12011 



20121 ). indicating it 



may be possible to detect such non-transiting objects by identifying these effects in their 
host star light curves. 

This paper presents the discovery of Kepler- 76b, the first hot Jupiter detected by the 
BEER algorithm that was subsequently confirmed by TRES and SOPHIE RV spectroscopy. 
It was identified by the BEER algorithm as a high-priority planetary candidate. Visual 
inspection of its light curve revealed a V-shaped primary transit and a minute secondary 
echpse, combined with beaming, ellipsoidal, and refiection amplitudes, consistent with a 
massive-planet com panion. We noticed later that this st ar was listed in the Kepler eclips- 
ing binary catalog (jPrsa et al.ll201ll : ISlawson et al.l |20 111 ). Based on this information spec- 
troscopic follow-up observations were initiated for this target, which in turn confirmed its 
planetary nature. 

Section 2 presents the BEER search and the initial analysis of the Kepler light curve. 
Section 3 provides the details and results of the spectroscopic observations. Section 4 de- 
scribes the details and results of the light curve transits and occultations analysis. Section 5 
presents the detection of evidence for superrotation in the light curve and Section 6 discusses 
the implications of, and conclusions from, the findings of this paper. 



2. The Photometric BEER search 

To identify candidates for low-mass companions we analyzed the Kepler raw light curves 
of the Q2 to QIO quarters, spanning 831 days. We visually identified 22 time segments that 
showed instrumental artifacts in the photometry, and ignored data points within those seg- 
ments, removing a total of 59.9 days of data from the light curves. We also corrected 
two systematic jumps at Kepler times (BJD — 2454833) of 200.32 and 246.19 days. For 
each light curve, outliers were then removed by 4a clipping and de trending was performed 



using a cosine-tran sform filter, adapted to unevenly spaced data (JMazeh &: Faiglerl 12010 
Mazeh et al.l |2012| ). resulting in a cleaned and detrended light curve. We then applied 



the BEER algorithm to 41,782 stars brighter than 13.7 mag, with Kepler Input Catalog 
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( iBrown et al.ll201ll ) radius smaller than ARq, calculating the Fast Fourier Transform (FFT) 
based power spectrum of the cleaned and detrended light curve of each star, interpolated 
over the gaps. Next, in order to avoid spurious peaks at long periods, we divided the full 
period range of each power spectrum into five sub-ranges: [0.3-1], [1-2], [2-5], [5-10], and 
[10-20] days, and identified the highest peak within each sub-range. For each of the five 
peaks we de rived the BEER amplit udes and the estimated mass and albedo of the candidate 
companion (JFaigler fc Mazehll2011l ). assuming the peak corresponds to either the orbital pe- 
riod (beaming and reflection) or half the orbital period (ellipsoidal). The BEE R amplitudes 



were c alculated by fitting the data with the modified BEER model suggested by lMazeh et al. 
(120 12l ). that uses a Lambertian reflection/emission function. 



We then selected 26 candidates with the highest signal to noise ratio for the ellipsoidal 
and beaming amplitudes, and with estimated secondary mass smaller than 60 Mjup and im- 
plied albedo smaller than 0.5. One of these candidates was Kepler- 76 {Kepler Input Catalog 
number 4570949), for which visual inspection revealed primary and secondary eclipses with 
depths of about 5 x 10"^ and 1 x 10~^, respectiv ely. This candida te was also identified as a 
member of the Kepler Eclipsing Binary catalog ( iPrsa et al.ll201ll ). Follow up spectroscopic 
observations confirmed the companion as a hot Jupiter. In forthcoming papers we will report 
on our observations of the 26 candidates, additional confirmation of a possible brown dwarf, 
and the false-positive rate of this sample. 

We report here the BEER analysis results for the light curve of Kepler-76 after mask- 
ing out the transit and occultation data points. We note that in our initial detection the 
BEER analysis was per formed on the unmasked data, but the use a robust-fit function 
( jHoUand fc Welschlll977l ) that identified the transit points as outliers resulted in similar mea- 
sured BEER amplitudes. Table [1] lists for Kepler-7 6 the coordinates and stellar properties 
from the Kepler Input Catalog (JBrown et al.l 120111 ) . revised effective temperature estimate 
from iPinsonneault et al.l ( 120121 ) , and the results of th e BEER analysis. Figu re [T] presents 
a short section of the 'cleaned' ( iMazeh fc Faiglerll2010l : iFaigler &: Mazehll201l[ ) photometric 
data of the host star. Figure [2] presents the FFT-based power spectrum, and Figure |3] shows 
the hght curve folded with the detected period. It is interesting to notice, by inspecting the 
cleaned light curves (Figured]) and the data and residuals r.m.s. (Tabled]), that the effects 
are significantly smaller than the light curve noise, to the point that the detected modula- 
tions almost can not be recognized by eye. However, deriving the BEER photometric power 
spectrum from data with time spans of hundreds of days, produces a prominent detectable 
peak at the orbital period (Figure |2]). 
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Table 1. Kepler- 76: Stellar properties and BEER results 



RA 


19:36:46. 11 


DEC 


39:37:08.4 


ftTp" [mag] 


13.3 


Teff" [K] 


6196 


TefF*" [K] 


6409 ± 95 


log 3"^ [dex] 


4.388 


[mlHY [dex] 


-0.033 


R." {R@\ 


1.12 


M.- [Mq] 


1.12 


h^ 


0.056 


BEER model: 




Period [days] 


1.5449 ±0.0007 


To - 2455000'' [BJD] 


737.49 ±0.19 


Ellipsoidal [ppm]'^ 


21.5 ±1.7 


Beaming [ppm]*^ 


15.6 ±2.2 


Reflection [ppm]"^ 


56.0 ±2.5 


r.m.s. cleaned data [ppm] 


133 


r.m.s. residuals [ppm] 


127 



Right ascension 

Declination 

Kepler band magnitude 

KIC Effective temperature 

Revised effective temperature 

Surface gravity 

Metallicity 

Primary radius 

Primary mass 

Average third light fraction 



With transit and occultation points masked out 

Orbital period 

Orbital zero phase ephemeris 

Ellipsoidal semi-amplitude 

Beaming semi-amplitude 

Lambcrtian reflection/emission semi-amplitude 

Root-mean-square of the cleaned and detrended data 

Root-mean-square of residuals of data from BEER model 



"from Kepler Input Catalog 
revised T^ff from iPinsonneault et al.l 1120121 ') 
'^calculated from Kepler Input Catalog log g and R 

''To is the time in which the companion is closest to the observer, assuming a circular orbit 
'^corrected for third light 
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Fig. 1. — The light curve of Kepler-76 for a selected time span of 37 days, after outlier 
removal and long-term detrending. Top: The untruncated light curve, showing the full 
depth of the transits. Bottom: The light curve with the core of the transit events truncated. 
Note the different scales of the two plots. (The transit missing at time 780.9 fell in a short 
gap in the raw data.) 
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Kepler-76 Power Spectrum 




Frequency [1/day] 



Fig. 2. — The FFT based power spectrum of the detection. The orbital period and half- 
orbital period peaks are marked by vertical dashed lines. For clarity, only the frequency 
range of 0-1.5 day~^ is plotted, since no significant peak was found for frequencies higher 
than 1.5 day^""^. 



x10" 



Kepler-76 




0.2 0.4 0.6 O.i 

-5 



.x10 



(0 

.-g 0- 

w 



»\s%^^^^^ tetjfilsSf li^ 



0.2 



0.4 



0.6 



Phase 



0.8 



Fig. 3. — The folded cleaned light curve binned into 100 bins. Phase zero is when the 
companion is closest to the observer, while phase 0.5 is when the primary is closest to the 
observer, assuming a circular orbit. The error bars represent la uncertainties of the median 
value of each bin, based on the scatter of data points within that bin. The line presents the 
BEER model. The residuals of the data from the model are plotted in the bottom panel. For 
clarity, the primary eclipse was removed. The secondary eclipse is clearly visible at phase 
0.5 of the plot. 
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3. Spectroscopic observations 

Spectroscopic observations of ttie candidate were obtained betwe en 29 May an d 6 Octo- 
ber 2012 with the Tillinghast Reflector Echelle Spectrograph (TRES; lFureszll2008[ ) mounted 
on the 1.5-ni Tillinghast Reflector at the Fred Lawrence Whipple Observatory operated by 
the Smithsonian Astrophysical Observatory (SAO) on Mount Hopkins in Southern Arizona, 
using the medium resolution flber at a spectral resolution of 44,000, covering a spectral 
range from 385 to 910 nm. Exposures of a Thorium-Argon hollow-cathode lamp immedi- 
ately before and after each exposure were used for wavelength calibration. The spectra were 
extracted and r ectifled to intensity vs . wavelength using standard procedures developed by 
Lars Buchhave (JBuchhave et al.ll2010l ). 



Additional spectroscopic obs ervations were obtained between 17 July and 1 Au gust 2012 
with the SOPHIE spectrograph JPerruchot et allboosl : ISouchv et all 120091 I2OI3I ) mounted 
on the 1.93-m telescope at Observatoire de Haute-Provence, France, using the High Efficiency 
mode {R ~ 39,000 at 550 nm) of the instrument. Spectra were extracted with the online 
standard pipeline. 



Following a method similar to the Stellar Parameter Classification method (SPC lBuchhave et al 



2OI2I ). the atmospheric parameters of Kepler- 76 were determined from the SOPHIE spec- 
tra by cross-correlating the observed spectral regio ns not affected by telluric lines against a 
library of synthetic spectra (JHauschildt et al.lll999f ). with varying values of effective temper- 
ature Teff, surface gravity log (7, metallicity [m/H] and rotational velocity vsini. For each of 
the observed spectra we derived the best set of parameters that yielded the highest correla- 
tion. This was done by fitting a second degree polynomial to the maximum correlation as 
function of each parameter around the synthetic spectrum that yielded the best correlation. 
The final parameter values for a star were taken as the mean of the parameter values derived 
for each observed spectrum of that star, weighted proportionally to the inverse of the scatter 
of the maximum around the fitted polynomial. 

The Phoenix library of synthetic spectra we used spans the following intervals in atmo- 
spheric parameters: 3000K < Teg < lOOOOK, -0.5 < \ogg < 5.5 (cgs), and -1.5 < [m/H] < 
+0.5. The spacing in T^s is lOOK for Teff < 7000K, and 200K elsewhere. The spacing in 
logg and [m/H] is 0.5 dex. The interval and spacing of vsini values in our algorithm are 
free parameters set by the user, since each s ynthet i c spec trum cho sen from the li b rary is 



convolved with a rotational profile G(v) (e.g., iGrayl (|2005[ ). p. 465; ISanterne et al.l (120 121 )) 



and a Gaussian representing the instrumental broadening of the lines, just before calculating 
cross- correlation function. 



The cross-correlation was performed using TODMOR (jZucker fc Mazelj|l994J : IZucker et al. 
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20031 . |2004| ) — a two-dimensional correlation algorithm, assuming the light contribution of 
the secondary is negligible. In TODMOR, the cross-correlation functions are calculated 
separately for each echelle order, and then corabined to a single cross- correlation function 
according the scheme proposed by IZucker et al.l (120031 ). The atmospheric parameters found 
this way are listed in Table [2l The relativel y large uncertainties result mainly from th e 
addition of possible systematic errors (see e.g., iBruntt et al.ll2010l l2012t iTorres et al.ll2012l ). 



The primary mass was estimated using the atmospheri c parameters derived fr om the 
spectra and a grid of Y^ stellar isochrones ( jYi et al.l l200ll : iDemarque et al.l l2004j ). This 
was done by taking into account all age and mass values that fall into the ellipsoid in 
the {Tes,logg, [Fe/H]) space defined by the atmospheric parameters and their errors. To 
illustrate the process Figure H] shows two sets of Y^ stellar isochrones of 0.2, 0.4, 1, 2, 4, 8, and 
10 Gyr — one for [Fe/H] = 0.05 (solid lines) and one for [Fe/H] = —0.27 (dashed lines). 
The ellipse defined by the estimated T^s and logg and their uncertainties is also shown. A 
lower limit of 0.2 Gyr on the stellar age was set to ignore possible pre-m ain sequence soluti ons. 
This procedure yielded a mass estimate of 1.20 ± 0.09Mq. Following iBasu et al.l (120121 ) we 
have conservatively doubled the mass errors to take into account possible uncertainties in 
stellar model parameters. 

Radial velocitie s were derived for t he TRES observations in two different ways, as de- 
scribed in detail by iFaigler et al.l ( 120121 ). First, we derived absolute velocities using cross- 
correlations of the observed spectra against the template from our library of synthetic spectra 
that yielded the best match (with Tgfj = 6000K, logg = 4.0 cgs, vsini = 12kms~^ and solar 
metallicity). The absolute velocity analysis used just the spectral order containing the Mgib 
triplet and was calibrated using lAU RV standard stars. With the goal of achieving better 
precision we also derived velocities using about two dozen spectral orders, correlating the 
individual observations against a template based on the strongest exposure. Thus the multi- 
order velocities are relative to the observation chosen as the template. They are reported in 
Table El 

For the SOPHIE observations, radial velocities were derived by computing the weighted 
cross-correlation function (CCF) of th e spectra with a numerical spectral mask of a G2V star 
( iBaranne et al.lll996l : lPepe et al.ll2002[ ). For the last five exposures, which wer e contaminated 



by sc attered moon light, we subtracted the sky using the fiber B spectrum ( ISanterne et al. 



20091 ). before deriving the radial velocities. Table |3] lists the radial- velocity measurements 
and their uncertainties. 
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Fig. 4. — Y^ stellar isochrones, from iDemarque et al.l ( 120041 ). of 0.2-10 Gyr for metallicities 
[Fe/H] = 0.05 (solid lines) and [Fe/H] = —0.27 (dashed lines). The estimated Tes and logg 
of Kepler-76 with their uncertainties are marked by a star and an ellipse. 



Table 2: Atmospheric parameters of Kepler-76 



Teff [K] 

log 5 [dex] 

[m/H] [dex] 
vsini [km s~^] 

M* [Mq] 



6300 ± 200 

4.2 ±0.3 

-0.1 ±0.2 

6.5 ±2 

1.2 ±0.2 



Table 3: Radial- velocity measurements 



Timo [BJD-2456000] RV [km b^'-] 



Instrument 



76.930366 

83.895818 

84.868623 

87.836880 

107.916759 

115.796005 

117.775394 

207.682884 

126.378842 

128.575439 

129.562543 

130.416149 

131.389038 

137.536798 

138.488304 

139.470933 

140.470516 

141.450747 



0.581 


0.069 


TRES 


0.161 


0.110 


TRES 


0.546 


0.072 


TRES 


0.607 


0.103 


TRES 


0.586 


0.098 


TRES 


0.727 


0.114 


TRES 


0.100 


0.082 


TRES 





0.069 


TRES 


4.999 


0.036 


SOPHIE 


5.597 


0.036 


SOPHIE 


5.081 


0.056 


SOPHIE 


5.560 


0.091 


SOPHIE 


5.341 


0.061 


SOPHIE 


5.196 


0.132 


SOPHIE 


5.194 


0.057 


SOPHIE 


5.615 


0.080 


SOPHIE 


4.992 


0.114 


SOPHIE 


5.171 


0.082 


SOPHIE 
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The first RV measurements of Kepler- 76 showed variability consistent with the pho- 
tometric orbital phase, so we continued observations in order to allow an orbital solution 
independent of the BEER analysis. The derived eccentricity of the solution was statistically 
indistinguishable from zero, so we reran the solution with eccentricity fixed to zero. Figure 
shows the follow-up RV measurements and the velocity curve for the orbital solution, folded 
with the period found, and the top section of Table H] lists the derived orbital elements for 
the independent RV solution. The center-of-mass velocities 7t and 7s for the independent 
RV sets from TRES and SOPHIE differ by 5.68 kms~^. This is because the TRES velocities 
are relative to the strongest observation, while the SOPHIE velocities are meant to be on 
an absolute scale. If the absolute TRES velocities derived using the Mg b order are used 
instead of the relative velocities, 7t = — 5.18kms~^, quite close to the SOPHIE value of 
7s = — 5.31kms~^. For the joint analysis reported below, the two independent velocity sets 
were shifted to a common zeropoint using the 7 velocities reported in Table HI 

Next, in order to obtain a combined solution from photometry and RV measurements, we 
reran the RV model using the photometric period and ephemeris, with their uncertainties, as 
priors. The bottom section of TablelUlists the orbital elements derived from this photometry- 
constrained RV solution, and the estimated minimum secondary mass, MpSini. 
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Fig. 5. — The RV measurements folded at the derived orbital period. In the top panel, the 
solid line presents the photometry-constrained orbital RV model and the horizontal-dashed 
line indicates the center-of-mass velocity. Circles denote SOPHIE RV points and squares 
denote TRES RV points. The residuals are plotted at the bottom panel. Note the different 
scales of the upper and lower panels. 
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Table 4. Orbital elements based on RV measurements 



N 


18 


Number of RV measurements 


span [days] 


130.8 


Time span of measurements 


Independent RV solution: 






To - 2456000 " [BJD] 


126.76 ±0.03 


Orbital cpliemeris 


P [days] 


1.5474 ± 0.0021 


Orbital period 


7t [km s-i] 


0.365 ± 0.042 


Zero point velocity of TRES 


7S [km s^-'-] 


-5.315 ±0.019 


Zero point velocity of SOPHIE 


Knv [km s-i] 


0.308 ± 0.020 


Radial velocity semi-amplitude 



RV solution constrained by the 
photometric period and ephemeris: 



To - 2456000 
P [days] 
7t [km s^-"-] 
7S [km s^l] 
Kbv [km s^l] 



[BJD] 



126.77 ±0.023 
1.5450 ±0.0005 

0.328 ± 0.032 
-5.320 ±0.020 

0.306 ± 0.020 



Orbital ephemeris 

Orbital period 

Zero point velocity of TRES 

Zero point velocity of SOPHIE 

Radial velocity semi-amplitude 



Mpsini [Mjup] 



1.96 ± 0.25 Minimum companion mass 



" To is the time in which the companion is closest to the observer 
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4. Photometric modeling of the hght curve 

For a more complete photometric analysis of this transiting hot Jupiter we used the 
Kepler light curves of the Q2 to Q13 quarters, spanning 1104 days. First we fitted the 
cleaned and detrended data with the BEER model while masking out data points in or 
around the transits and occultations. The fitted amplitudes, after correction for a third 
light using the KIC estimate, are listed in Table El We then subtracted the BEER model 
from the data and analyzed the data points in and around the transits and occultations. 
For that we ran a Markov chain Monte Car lo (MCMC) ana l ysis, w hile fitting the transit 



data points using a long-cadence integrated iMandel fc Agoll (120021 ) model with quadratic 



limb darkening, assuming a circular orbit. The model limb darke ning coefficients could not 



be constrained, so we kept them fixed at values interpolated from lClaret fc BloemenI (120111 ) 
using the stellar parameters derived from spectroscopy. We then fitted the occultation data 
keeping the geometric parameters derived from the transit fixed, and assuming a linear limb 
darkening coefficient of 0.5 for the planet, while looking for the occultation depth that best 
fits the data. 



Figure m presents the cleaned and detrended data points and the best-fit iMandel fc Agol 



(120021 ) model combined with the BEER model, both folded at the orbital period. Table |5llists 
the priors and the MCMC medians and la uncertainties of the posterior parameters. We note 
here that the orbital period P and time of primary transit Tq listed in Table El were derived 
from the transit data points, thus yielding high-accuracy estimates. The same parameters, 
listed in Table [H were derived from the out-of-transit data, thus yielded much lower- accuracy 
estimates. The To values listed in both tables were not corrected for the Kepler timing error 
http://archive.stsci.edu/kepler/timing_error.html), which should be taken in account when 



comparing them to non-Kepler observations. 
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Fig. 6. — Kepler-76 cleaned and detrended data points and best-fit model, both folded at 
the orbital period. In the top panel of each plot, the solid line presents the best-fit model, 
and the dots present the data points. The residuals are plotted in the bottom panel. Note 
the different scales of the upper and lower panel of each plot. 
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Table 5. Photometric light curve transit and occupation model 



Parameter Median and la uncertainty 



Priors: 

M* [Mq] 

a* 

6* 

dp 

6p 

/3 



1.2 ± 0.2 Primary mass 

0.313 Star linear limb darkening coefficient 

0.304 Star quadratic limb darkening coefficient 

0.5 Planet linear limb darkening coefficient 

Planet quadratic limb darkening coefficient 

0.056 Average third light fraction (From KIC) 



Posteriors: 

To [BJD] 
Period [days] 
i [deg] 
R*/a 
Rp/a 
d2 [ppm] 



2454966.54811 ±0.00007 

1.54492875 ±0.00000027 

78.0 ±0.2 

0.221 ± 0.003 

0.0214 ± 0.0008 

98.9 ±7.1 



Time of primary transit 
Orbital period 
Orbital inclination 
Fractional primary radius 
Fractional planet radius 
Occultation depth 



p*/pe 
Rt [Rq] 



Rp [R. 



JupJ 



0.522 ± 0.019 Primary density relative to the Sun density 

1.32 ± 0.08 Primary radius 

1.25 ± 0.08 Planet radius 

0.944 ±0.011 Impact parameter 
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5. Inflated beaming amplitude and planet equatorial super-rotating jet 

The spectroscopic RV observations and the hght curve transit and occultation analysis 
yielded independent orbital solutions with ephemeris and period nicely consistent with the 
BEER ephemeris and period. To compare the measured RV amplitude with the beaming- 
based predicted RV amplitude, one needs to evaluate the abeam factor, which corrects for 
the Doppler shift of the s t ellar spectrum relative to the observed band of the telescope 
jRvbicki fc LightmanlEoTol : IPaider fc Mazehlboilh . 



To estimate the abeam value, we used spectra from the library of ICastelli fc Kurucz 



(120041 ) models close to the estimated temperature, metallicity and gravity of the primary 
star, numerically shifting them relative to the Kepler response funct i on, while taking into 



account the photo n counting nature of Kepler (JLoeb &: Gaudil l2003l : iBloemen et al.l 12011 



Faigler et al 



20121 ) ■ For c larity we note that by definition ab 



4 



where a is the 

power-law index used bv iLoeb fc Gaudil (2003) a r id jB ) is the photon weighted bandpass- 
integrated beaming factor used by lBloemen et al.l (J201ll ). The result of this calculation gave 
abeam = 0.92 ± 0.04, resulting in an RV semi-amplitude of -ft'beam = l-H ± 0.17kms~^. The 
RV semi-amplitude predicted from the beaming effect was 3.5 times larger than the measured 
amplitude, with a difference significance of about 4.5o" between the two. 

A possible explanation for this inflated photometric beaming amplitude might be a phase 
shift of the reflection signal, due to the superrotation phenomenon, which involves eastward 
advection of gas by an equatorial s uper- rotating jet within the atmosphere of a co-rotating 
companion. IShowman fc GuillotI (120021 ) predicted through a 3D atmospheric circulation 



model that tidally-locked, short period planets develop a fast eastward, or superrotating, jet 
stream that extends from the equator to latitudes of typically 20° — 60°. They showed that 
in some cases (depending on the imposed stellar heating and other factors) this jet causes an 
eastward displacement of the hottest regions by 10° — 60° longitude from the substellar point, 
resulting in a pha se shift of the thermal emiss ion phase curve of the planet. This prediction 
was confirmed by iKnutson et al.l (120071 . |2009| ) through Spitzer infrared observations of HD 
189733, which indicated a phase shift of 16° ± 6° in the 8 /im band and 20° — 30° in the 
24 /im band. 

In general, what we call a reflection modulation is actually the light scattered off the 
planet in combination with radiation absorbed and later thermally re-emitted at different 
wavelengths. The two processes are controlled by the Bond albedo, < A^ < 1, and the 
day-night heat redistribution efficiency, < e < 1, which can be constrained only if observa- 
tions of the phase mo dulation or the secondary eclipse are available in different wavelengths 
( ICowan fc Agolll201ll ). This makes it impossible to distinguish between r eflected and re- 
radiated photons from the single band Kepler light curve we have in hand. ICowan fc Agol 
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( 120111 ) discuss HAT-P-7 as an example, and show that its Kepler light curve can be explained 
as mostly reflected light at one limit, to mostly thermal emission at the other limit, with an 
entire range of models between them being consistent with the light curve. This is important 
for the current discussion, as we expect superrotation to shift only the thermal re-emission, 
while leaving the scattered light component unshifted. 

To estimate the maximum frac tion of the reflection a mplitude originating from thermal 



re-emission in our case, we follow I Cowan fc Agoll (J201ll ) and estimate the no albedo, no 



redistribution, effective day side temperature T^=o ~ 2670/^, which translates in the Kepler 
band to a maximum reflection amplitude Aj-ef ^ 37 ppm. This means that the measured 
amplitude of ~ 50 ppm can be explained mostly by thermal re-emission. The actual fraction 
of thermal emission in this case is probably smaller, but this calculation illustrates that 
the fraction of thermal emission in the visual Kepler light curve phase modulation may be 
significant, making it a worthy effort to look for a superrotation phase shift in the light curve. 

We suggest here, that if such a phase shift is present in the Kepler light curve, it will 
show up in our phase curve model mainly as an inflated beaming amplitude. To illustrate 
that we consider a simple superrotation model consisting of a phase-shifted geometric reflec- 
tion/emission combined with a beaming modulation, 

MsR = -^ref COS(0 + 6sr) + ^beam slu = -A^ei COS 6sR COS (j) + (Abeam + ^ref slu 6sr ) slu (j) , 

(1) 

where Aref is the reflection/emission semi-amplitude, Abeam is the beaming semi-amplitude, 
is the orbital phase relative to mid-transit, and 6sr is the superrotation phase shift angle. 
This model suggests that if a phase shift is present and the reflection amplitude is larger than, 
or of the order of, the beaming amplitude, the underlined term in Eq. 1 may add substantially 
to the amplitude of the sin modulation, mimicking an inflated beaming effect. 

To test our conjecture that the beaming/ellipsoidal inconsistency is a result of a su- 
perrotation phase shift of the reflection/emission phase modulation, we fitted the data 
using the derived syste m parameters (Tables |5] and [2]) and the BEER effects equations 



( JFaigler fc Mazehll201l[ ). while looking for the planetary mass, geometric albedo, and phase 
shift of the Lambertian phase function that minimized the x^ of the fit. Adding the phase 
shift parameter to the model resulted in a decrease of the x^ value by 90, relative to the no- 
phase-shift model, indicating a substantially better agreement of the data with a model that 
combines beaming, ellipsoidal and a phase shifted Lambertian reflection. An F-test shows 
that fitting the data while allowing for a phase shift, as opposed to the no-phase- shift null 
model, yields a better fit with a confidence level better than 9a. Table |6] lists the amplitudes 
derived by the BEER analysis, the planetary mass derived directly from the beaming versus 
the ellipsoidal amplitudes, and the spectroscopic RV derived planetary mass. The table then 
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lists the planetary mass, phase-shift angle and geometric albedo resulting from the superro- 
tation model. The superro tation phase shift est i mate is small and well within the theoretical 
limit of 60° predicted by IShowman &: GuillotI (120021 ) . In addition, the derived planetary 
mass estimate is well within the la range of the RV measured planetary mass, indicating 
that assuming superrotation resolves the inconsistency and provides a good estimate for the 
planetary mass, derived solely from the Kepler photometry, given a good stellar model. 
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Table 6. BEER amplitudes, planetary mass, and superrotation phase shift angle 



BEER model: 
Ellipsoidal [ppm]" 
Beaming [ppm]" 
Reflection [ppm]" 
a2s [ppm]" 



21.1 it 1.7 Ellipsoidal semi-amplitude 

13.5 ± 2.0 Beaming semi-amplitude 

50.4 ± 2.0 Reflection semi-amplitude 

—2.7 ± 1.0 sin2ij> coefficient 



Mp,bcam [Mjup] 
Mp,ellip [A^Jup] 
Mp,RV [Mjup] 



7.2 ± 1.4 Planetary mass derived from beaming 
2.1 ± 0.4 Planetary mass derived from ellipsoidal 
2.00 ± 0.26 Planetary mass measured by spectroscopic RV 



Superrotation phase-shift model solution: 
Mp,SR [Mjup] 
5sR [deg] 
A, 



N 

^null 



2.1 ±0.4 

10.3 ±2.0 

0.23 ±0.02 

31468 

31465 

31555 



Planetary mass 

Superrotation phase shift angle 

Geometric albedo 

Number of data points 

X^ of a model allowing for a phase shift 

X^ of a zero phase shift model (null model) 



"corrected for third light 
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Discussion 



This paper presents a new hot- Jupiter companion, Kepler- 76b, initially identified by 
the BEER algorithm, and later confirmed by spectroscopic observations. The BEER de- 
tection was based on the photometrically measured amplitudes of the BEaming, Ellipsoidal 

and Refiection effects, that were consistent with a planeta ry companio n. This is just the 

third confirmed plan et in the Kepler Held, after HAT-P-7b ( IWelsh et al.ll2010l ) and TrES-2b 
( JBarclay et al.l 120121 ). that its host light curve exhibits the three phase curve effects, and is 
the faintest of the three stars. It is also one of a few confirmed grazing exoplanets, showing 
a V-shaped transit and a partial occultation. 

We have identified an inconsistency between the beaming amplitude and the spectro- 
scopically measured radial velocity. Similar inconsistencies between the planetary mass de- 
rived from the beaming amplitudes and the mass derived from the ellipsoidal amplitude were 
noticed previously by several authors for KOI-13 (JMazeh et al.ll2012l : IShporer et al.ll201l[ ) and 
TrES-2 f lBarclay et al.ll2012[ ). We suggest here that these inconsistencies can be explained 
by a phase shift of the p lanetary thermal modulatio n due to the equatoria l superrotation 



phenomena predicted by IShowman fc GuillotI ( 12002| ) and later observed by iKnutson et al. 



( 120071 . I2OO9I ) in the infrared for HD 189733. In such cases we should be able to measure the 
superrotation phase-shift angle from the visual band Kepler light curve of the system. As 
we do not expect scattered light to exhibit such a phase shift, visual band detection of the 
superrotation phase shift may yield a constraint on the ratio of scattered light to thermally 
re-emitted light from the planet. 

Finally, we wish to briefiy comment on the sensitivity of the BEER algorithm. The 
detection presented here exhibits the lowest-mass companion identified so far by the algo- 
rithm, indicating its possible current detection limit. As we require that a BEER candidate 
must show statistically significant beaming and ellipsoidal effects to be considered valid, 
we choose here the minimum of the semi-amplitudes of the two effects as the BEER de- 
tectability parameter of a planet. To estimate our ability to detect more p lanets and brown 
dwarf s. Figure [7] presents the calculated value of this parameter, using the iFaigler fc Mazeh 
(I2OIII ) equations, for known exoplanets of mass higher than 0.5Mjup and period shorter than 



30 days, as of January 2013 (]http://exoplanet.eu/ ), together with the measured value for 
Kepler-76b. The figure shows that there are six transiting and one RV detected planets with 
calculated amplitudes higher than that of Kepler- 76b, suggesting that these systems could 
have been detected by the BEER algorithm, if their stellar and instrumental noise were 
similar to that of Kepler- 76. It is also apparent that further enhancement of the algorithm 
sensitivity could significantly increase the number of potentially detectable planets. Given 
the fact that BEER can detect similar non-transiting objects, we expect to find more objects 
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once we improve our detection threshold. 
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Fig. 7. — The minimum of the predicted beaming and elhpsoidal semi-amphtudes, as a 
function of the secondary mass, for known exoplanets of mass higher than O.SMjup and 
period shorter than 30 days. Blue dots mark exoplanets detected by RV, red circles mark 
exoplanets detected by transit, and the red star marks the measured amplitude of Kepler-76b. 
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